Bioinformatics tools for analysing viral genomic data

The field of viral genomics and bioinformatics is experiencing a strong resurgence due to high-throughput sequencing (HTS) technology, which enables the rapid and cost-effective sequencing and subsequent assembly of large numbers of viral genomes. In addition, the unprecedented power of HTS technologies has enabled the analysis of intra-host viral diversity and quasispecies dynamics in relation to important biological questions on viral transmission, vaccine resistance and host jumping. HTS also enables the rapid identification of both known and potentially new viruses from field and clinical samples, thus adding new tools to the fields of viral discovery and metagenomics. Bioinformatics has been central to the rise of HTS applications because new algorithms and software tools are continually needed to process and analyse the large, complex datasets generated in this rapidly evolving area. In this paper, the authors give a brief overview of the main bioinformatics tools available for viral genomic research, with a particular emphasis on HTS technologies and their main applications. They summarise the major steps in various HTS analyses, starting with quality control of raw reads and encompassing activities ranging from consensus and de novo genome assembly to variant calling and metagenomics, as well as RNA sequencing

Analysis of Paramyxovirus Transcription and Replication by High-Throughput Sequencing.

We have developed a high-throughput sequencing (HTS) workflow for investigating paramyxovirus transcription and replication. We show that sequencing of oligo(dT)-selected polyadenylated mRNAs, without considering the orientation of the RNAs from which they had been generated, cannot accurately be used to analyze the abundance of viral mRNAs because genomic RNA copurifies with the viral mRNAs. The best method is directional sequencing of infected cell RNA that has physically been depleted of ribosomal and mitochondrial RNA followed by bioinformatic steps to differentiate data originating from genomes from viral mRNAs and antigenomes. This approach has the advantage that the abundance of viral mRNA (and antigenomes) and genomes can be analyzed and quantified from the same data. We investigated the kinetics of viral transcription and replication during infection of A549 cells with parainfluenza virus type 2 (PIV2), PIV3, PIV5, or mumps virus and determined the abundances of individual viral mRNAs and readthrough mRNAs. We found that the mRNA abundance gradients differed significantly between all four viruses but that for each virus the pattern remained relatively stable throughout infection. We suggest that rapid degradation of non-poly(A) mRNAs may be primarily responsible for the shape of the mRNA abundance gradient in parainfluenza virus 3, whereas a combination of this factor and disengagement of RNA polymerase at intergenic sequences, particularly those at the NP:P and P:M gene boundaries, may be responsible in the other viruses.IMPORTANCE High-throughput sequencing (HTS) of virus-infected cells can be used to study in great detail the patterns of virus transcription and replication. For paramyxoviruses, and by analogy for all other negative-strand RNA viruses, we show that directional sequencing must be used to distinguish between genomic RNA and mRNA/antigenomic RNA because significant amounts of genomic RNA copurify with poly(A)-selected mRNA. We found that the best method is directional sequencing of total cell RNA, after the physical removal of rRNA (and mitochondrial RNA), because quantitative information on the abundance of both genomic RNA and mRNA/antigenomes can be simultaneously derived. Using this approach, we revealed new details of the kinetics of virus transcription and replication for parainfluenza virus (PIV) type 2, PIV3, PIV5, and mumps virus, as well as on the relative abundance of the individual viral mRNAs

Genome sequences of five African swine fever virus genotype IX isolates from domestic pigs in Uganda

Complete genome sequences of five African swine fever virus isolates were determined directly from clinical material obtained from domestic pigs in Uganda. Four sequences were essentially identical to each other, and all were closely related to the only known genome sequence of p72 genotype IX

Genome sequences of five African swine fever virus genotype IX isolates from domestic pigs in Uganda

Complete genome sequences of five African swine fever virus isolates were determined directly from clinical material obtained from domestic pigs in Uganda. Four sequences were essentially identical to each other, and all were closely related to the only known genome sequence of p72 genotype IX

High Throughput siRNA Screening Identifies Phosphatidylinositol 3-kinase Class II Alpha as Important for Production of Human Cytomegalovirus Virions.

High throughput siRNA screening is a useful methodology to identify cellular factors required for virus replication. Here we utilized a high throughput siRNA screen based on detection of a viral antigen by microscopy to interrogate cellular protein kinases and phosphatases for their importance during human cytomegalovirus (HCMV) replication, and identified the Class II Phosphatidylinositol 3-kinase PI3K-C2A as being involved in HCMV replication. Confirming this observation, infected cells treated with either pooled or individual siRNAs targeting PI3K-C2A mRNA produced approximately 10-fold less infectious virus compared to controls. Western blotting and quantitative PCR analysis of infected cells treated with siRNAs indicated that depletion of PI3K-C2A slightly reduced accumulation of late, but not immediate-early or early, viral antigens and had no appreciable effect of viral DNA synthesis. Analysis of siRNA treated cells by electron microscopy and western blotting indicated that PI3K-C2A was not required for production of viral capsids, but did lead to increased numbers of enveloped capsids in the cytoplasm that had undergone secondary envelopment and reduction of viral particles exiting the cell. Therefore, PI3K-C2A is a factor important for HCMV replication and has a role in production of HCMV virions. IMPORTANCE: There is limited information about the cellular factors required for human cytomegalovirus (HCMV) replication. Therefore, to identify proteins involved in HCMV replication we developed a methodology to conduct a high throughput siRNA screen in HCMV infected cells. From our screening data we focused our studies on the top "hit" from our screen, the lipid kinase phosphatidylinositol 3-kinase Class II Alpha (PI3K-C2A), as its role in HCMV replication was unknown. Interestingly, we found that PI3K-C2A is important for the production of HCMV virions and is involved in virion production after secondary envelopment of viral capsids, the encapsidation of HCMV capsids by a lipid bilayer that occurs before virions exit the cell

Development of a reverse genetics system for Toscana virus (lineage A)

Toscana virus (TOSV) is a Phlebovirus in the Phenuiviridae family, order Bunyavirales, found in the countries surrounding the Mediterranean. TOSV is an important cause of seasonal acute meningitis and encephalitis within its range. Here, we determined the full sequence of the TOSV strain 1500590, a lineage A virus obtained from an infected patient (Marseille, 2007) and used this in combination with other sequence information to construct functional cDNA plasmids encoding the viral L, M, and S antigenomic sequences under the control of the T7 RNA promoter to recover recombinant viruses. Importantly, resequencing identified two single nucleotide changes to a TOSV reference genome, which, when corrected, restored functionality to the polymerase L and made it possible to recover infectious recombinant TOSV (rTOSV) from cDNA, as well as establish a minigenome system. Using reverse genetics, we produced an NSs-deletant rTOSV and also obtained viruses expressing reporter genes instead of NSs. The availability of such a system assists investigating questions that require genetic manipulation of the viral genome, such as investigations into replication and tropism, and beyond these fundamental aspects, also the development of novel vaccine design strategies

Human cytomegalovirus genomes sequenced directly from clinical material: variation, multiple-strain infection, recombination and gene loss

The genomic characteristics of human cytomegalovirus (HCMV) strains sequenced directly from clinical pathology samples were investigated, focusing on variation, multiple-strain infection, recombination, and gene loss. A total of 207 datasets generated in this and previous studies using target enrichment and high-throughput sequencing were analyzed, in the process enabling the determination of genome sequences for 91 strains. Key findings were that (i) it is important to monitor the quality of sequencing libraries in investigating variation; (ii) many recombinant strains have been transmitted during HCMV evolution, and some have apparently survived for thousands of years without further recombination; (iii) mutants with nonfunctional genes (pseudogenes) have been circulating and recombining for long periods and can cause congenital infection and resulting clinical sequelae; and (iv) intrahost variation in single-strain infections is much less than that in multiple-strain infections. Future population-based studies are likely to continue illuminating the evolution, epidemiology, and pathogenesis of HCMV

Distemper, extinction, and vaccination of the Amur tiger

Canine distemper virus (CDV) has recently emerged as an extinction threat for the endangered Amur tiger (Panthera tigris altaica). CDV is vaccine-preventable, and control strategies could require vaccination of domestic dogs and/or wildlife populations. However, vaccination of endangered wildlife remains controversial, which has led to a focus on interventions in domestic dogs, often assumed to be the source of infection. Effective decision making requires an understanding of the true reservoir dynamics, which poses substantial challenges in remote areas with diverse host communities. We carried out serological, demographic, and phylogenetic studies of dog and wildlife populations in the Russian Far East to show that a number of wildlife species are more important than dogs, both in maintaining CDV and as sources of infection for tigers. Critically, therefore, because CDV circulates among multiple wildlife sources, dog vaccination alone would not be effective at protecting tigers. We show, however, that low-coverage vaccination of tigers themselves is feasible and would produce substantive reductions in extinction risks. Vaccination of endangered wildlife provides a valuable component of conservation strategies for endangered species

Highly Diverse Hepatitis C Strains Detected in Sub‐Saharan Africa Have Unknown Susceptibility to Direct‐Acting Antiviral Treatments

The global plan to eradicate hepatitis C virus (HCV) led by the World Health Organization outlines the use of highly effective direct‐acting antiviral drugs (DAAs) to achieve elimination by 2030. Identifying individuals with active disease and investigation of the breadth of diversity of the virus in sub‐Saharan Africa (SSA) is essential as genotypes in this region (where very few clinical trials have been carried out) are distinct from those found in other parts of the world. We undertook a population‐based, nested case‐control study in Uganda and obtained additional samples from the Democratic Republic of Congo (DRC) to estimate the prevalence of HCV, assess strategies for disease detection using serological and molecular techniques, and characterize genetic diversity of the virus. Using next‐generation and Sanger sequencing, we aimed to identify strains circulating in East and Central Africa. A total of 7,751 Ugandan patients were initially screened for HCV, and 20 PCR‐positive samples were obtained for sequencing. Serological assays were found to vary significantly in specificity for HCV. HCV strains detected in Uganda included genotype (g) 4k, g4p, g4q, and g4s and a newly identified unassigned g7 HCV strain. Two additional unassigned g7 strains were identified in patients originating from DRC (one partial and one full open reading frame sequence). These g4 and g7 strains contain nonstructural (ns) protein 3 and 5A polymorphisms associated with resistance to DAAs in other genotypes. Clinical studies are therefore indicated to investigate treatment response in infected patients. Conclusion: Although HCV prevalence and genotypes have been well characterized in patients in well‐resourced countries, clinical trials are urgently required in SSA, where highly diverse g4 and g7 strains circulate.Supported by the Medical Research Council (MRC) (MC_UU_12014/1) and Wellcome Trust (102789/Z/13/A) (to E.T.). M.S. is funded by the Wellcome Trust Sanger Institute (WT098051), the National Institute for Health Research Cambridge Biomedical Research Centre, the African Partnership for Chronic Disease Research (MRC UK partnership grant number MR/K013491/1), and the UK MRC (G0901213‐92157, G0801566). P.K. is funded by the UK MRC and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement. J.S. is funded by the MRC Confidence in Concept award to the University of Glasgow (MC PC 16045). G.M. is a Gates Cambridge Scholar supported by the Gates Cambridge Trust